spectroscopic parameters for silacyclopropynylidene, SiC2, from extensive astronomical 
observations toward CW Leo (IRC +10216) with the Herschel satellite^''^'^ 
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Abstract 



A molecular line survey has been carried out toward the carbon-rich asymptotic giant branch star CW Leo employing the HIFI 
instrument on board of the Herschel satellite. Numerous features from 480 GHz to beyond 1 100 GHz could be assigned unambigu- 
,— I ously to the fairly floppy SiCa molecule. However, predictions from laboratory data exhibited large deviations from the observed 
frequencies even after some lower frequency data from this survey were incorporated into a fit. Therefore, we present a combined 
fit of all available laboratory data together with data from radio-astronomical observations. 
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1. Introduction 

Silacyclopropynylidene, SiCa, somewhat better known as sil- 
icon dicarbide, is a fascinating molecule for spectroscopists, 
structural and quantum chemists as well as astronomers. In 
1926, uncataloged bands near 500 nm were discovered in the 
spectra of several carbon-rich asymptotic giant branch (AGB) 
stars m 01 ■ These are late-type stars which produce ele- 
ments heavier than helium and which eject large quantities of 
'gaseous material as well as dust which form a circumstellar en- 
velope (CSE). Thirty years later, laboratory spectroscopy estab- 
lished that the molecule SiCa is the carrier of these bands ||3|l. 
It was assumed that the molecule has a linear SiCC struc- 
ture (silapropadienediylidene) in analogy with the isoelectronic 
propadienediylidene, C3. Although some later studies cast 
doubt on the linear structure of the molecule, it took almost 
another 30 years until the analysis of the rotational structure of 
the electronic origin band unmistakably determined the struc- 
ture as silacyclopropynylidene [4]. In the course of their analy- 
sis, the authors instigated quantum chemical calculations which 
provided evidence that the cyclic isomer of SiC2 may be lower 
in energy than the linear form Js]] . 



We dedicate this work to tlie memory of Gisbert Winnewisser, a pioneer of 
Herschel and of astrocliemistry. 

** Herschel is an ESA space observatory with science instruments provided 
by European-led Principal Investigator consortia and with important participa- 
tion from NASA. 
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A plethora of quantum chemical calculations on various 
properties of SiC2 have been published later, yielding energy 
differences between the linear and the cyclic form which de- 
pended strongly on the level of the calculation and the size of 
the basis set. A high level ab initio calculation concluded that 
the cyclic isomer of SiCa is the only minimum on the potential 
energy surface and that the linear transition state is 24.3 kJ/mol 
higher in energy |6]. However, the authors attached a caveat to 
this value as an anharmonic force field calculation provided a 
much too small value for the vibrational energy of LI3 = 1 and 
much too large anharmonicity constants. Therefore, some of 
the authors revisited the problem of the energy difference be- 
tween the two SiC2 structures, the last time in 2003 when very 
high level calculations combined with very large basis sets, ba- 
sis set extrapolation to infinite size as well as additional correc- 
tions yielded a value of 26.5 kJ/mol [7]. 

The SiC2 molecular parameters obtained in Ref. |4i] laid 
the foundation for progress in laboratory spectroscopy. The 
7=1-0 rotational transition frequencies of the three iso- 
topologs Sicili ^'^SiCa, and ^"SiCa, as well as the permanent 
electric dipole moment were measured using Fourier transform 
microwave spectroscopy [8]. Subsequently, 34 additional tran- 
sition frequencies were measured for the main isotopic species 
between 93 and 370 GHz |@]. Even though a comparatively 
large number of 15 spectroscopic parameters of a standard 
Watson-type Hamiltonian in the A-reduction, all parameters up 
to sixth order, were employed in the fit, the transition frequen- 
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cies were reproduced on average to only four times the experi- 
mental uncertainties. 

Similarly large bodies of laboratory transition frequencies 
were obtained for SiC'^^C between 339 and 405 MHz [10] and, 
only very recently, for ^^SiCa and ^"SiCa between 140 and 
360 GHz fJ3i. 

Higher excited vibrational states have been studied for the 
main isotopolog. Rotational transitions in its low-lying vj, - I 
vibrational state between 186 and 399 GHz [T^ and in 1^3 = 
1 and 2 were obtained between 140 and 400 GHz (l^. The 
symmetry of - I is b2 and can be viewed as an asymmetric 
bending state which facilitates internal rotation of the C2 unit 
with respect to the Si atom. Its vibrational energy has been 
determined as 196.37 cm ' from an investigation into the laser- 
induced and the dispersed fluorescence of a jet-cooled sample 
ofSiCz 114J. 

Several attempts have been made to model rotational and 
sometimes also rovibrational data to a varying extent and ac- 
curacy. A semirigid bender (SRB) Hamiltonian was em- 
ployed lfl6ll to rep roduce - and 1 rotational transition 
frequencies [SO as well as rovibrational data from their pre- 
vious ifTill and present work |16]. The analysis suggested that 
the Unear configuration is not a local minimum, and the cyclic 
form is 22.5 kJ/mol lower than the linear form; the estimated 
uncertainty was 2.4 kJ/mol. They found the energy difference 
to be in good agreement with results from their own ab initio 
calculations; their highest level value being 21.8 kJ/mol. The 
reproduction of the vibrational data was reasonable, that of the 
rotational data was only qualitative and thus of no use for radio- 
astronomical purposes. 

The V3 - and 1 state rotational data [l2l were also mod- 
eled with a dedicated internal rotation Hamiltonian ifisll . A 
greatly improved reproduction, within 1.5 times the experimen- 
tal uncertainties, was achieved for the ground vibrational state, 
albeit at the expense of 16 spectroscopic parameters compared 
to 15 previously |9]. The barrier to linearity was derived as 
~54 kJ/mol, more than twice the value from ab initio calcula- 
tions of that time. 

A reasonably successful result has been obtained by em- 
ploying a conventional Watson-type Hamiltonian in the 5- 
reduction iHtIi . The reproduction of the ground state rotational 
data lUl^] was converged after varying 17 spectroscopic and 
keeping one fixed. The data were reproduced on average as 
well as in the fit which employed an internal rotation Hamilto- 
nian 11511 ; the number of varied parameters, however, was larger 
by still another one. After scaling of the parameters, the tran- 
sition frequencies of ^^SiCa and ^*'SiC2 available at that time, 
mostly from astronomical observation could be repro- 

duced well after releasing only 5 spectroscopic parameters 1 171 ; 
13 parameter were released in the fit for SiC'^C because of the 
large amount of accurate laboratory data fl^ . 

The unambiguous assignment of a spectroscopic feature ob- 
served in space relies on reliable predictions which are usu- 
ally based on laboratory data which, for the most part, have 
been obtained in approximately the same frequency domain. 
The predictions may also be based in part on data from astro- 
nomical observations, recent examples include H'-'CO^ 



DCO+ yj I20I], DNC and HN'^C |l20[], C2H UM, and 
C'^CH 02211 . Sometimes, identifications are even possible in 
the absence of laboratory data, as demonstrated recently by the 
detection of C5N in the CSE of CW Leo |[23J. Frequencies 
have to be predicted with accuracies better than around one 
tenth of the line width to permit extraction of dynamical in- 
formation. Some type of intensity information, such as that at a 
certain temperature, the line strength, or the Einstein A-value is 
needed as is additional auxiliary information such as quantum 
numbers, lower or upper state energies etc. Pickett's spcat and 
SPFIT programs ll24ll have been developed for that purpose and 
have evolved over time, see e. g. Ref. |25|. These programs are 
routinely used in the Cologne Database for Molecular Spec 
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troscop>t3, CDMS |2^ 27], to provide in its catalog sectiord 
predictions of (mostly) rotational spectra of molecules which 
may be found in various environments in space. 

The spectroscopic parameters from the first rotational analy- 
sis of the electronic spectrum of SiCa fT] were accurate enough 
to identify nine unassigned emission features, previously ob- 
served between 93 and 171 GHz toward the carbon-rich AGB 
star CW Leo, also known as IRC -1-10216, as belonging to SiCa 
and improved the SiCa structural parameters tZS\\ . Using the 
structural parameters from that work, the three J = 4-3, 
AKa - transitions of ^^SiCa and ^"SiCa were detected to- 
ward the same source very soon thereafter [29] ; two of these 
transitions were marginally detected for SiC'^C. Extensive sets 
of transition frequencies of ^^SiCa, ^"SiCa, and SiC'-'C were 
obtained from radio-astronomical observation between 90 and 
241 GHz together with laboratory rest frequencies for the latter 
isotopolog fl(J\. 

Radio lines of SiCi have also been detected in the CSEs 
of other C-rich AGB stars, e. g. toward II Lup, which is also 
known as IRAS 15194-5115 [^]. However, SiCa features 
are particularly strong toward CW Leo, which, to a large ex- 
tent, is due to its proximity to our Solar system. CW Leo has 
been studied extensively, and many molecular species, such as 
CN" [31] and FeCN [32J, have been detected toward its CSE 
exclusively or for the first time, including many Si-containing 
ones, such as SiCN and SiNC (3^. The source is a sub- 
ject of observations in several key projects carried out with 
the recently launched Herschel satellite 1341] . One of these 
projects is a molecular line survey carried out with the Het- 
erodyne Instrument for the Far-Infrared (HIFI) [35,] . This high- 
resolution instrument covers, in several bands, the 480-1250 
and 1410-1910GHz regions. A preliminary analysis ofalower 
frequency region (554-637 GHz) revealed that a rather small 
number of molecules account for a large fraction of the emis- 
sion features f^^. These molecules are CS, SiO, SiS, HCN, 
and, in particular, SiCa. Predictions of the rotational spectrum 
based on Ref. lll7[l turned out to be very good for transitions 
with higher values of Ka, but showed increasing deviations of 
up to 10 MHz for transitions with decreasing Ka, contrary to 
common expectations. A combined fit of these transitions fre- 
quencies together with those from laboratory spectra [HI ^ re- 
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quired only one additional parameter 113611 to reproduce the as- 
tronomical data within uncertainties, and the laboratory data 
about as well as before [17]. Predictions of the SiC2 rotational 
spectrum based on these results are currently available as ver- 
sion 2 in the CDMS0. Since they were still not appropriate to 
predict the observed emission features satisfactorily to high fre- 
quencies (beyond 1100 GHz), we present here a combined fit 
using these as well as laboratory data, supplemented with data 
from additional astronomical observations. 



2. Observed spectrum 

Silacyclopropynylidene, SiCa, is a triangular, fairly asym- 
metric rotor with k ^ {2B - A - C)/(A - C) = -0.7117. Its 
dipole moment of 2.393 (6) D Jst] is along the a-axis. Only 
even Kg exist because of the two equivalent C nuclei with zero 
spin. -branch transitions (A7 - J' - J" - +Y) with AKa - 
are the strongest ones, and most of the transitions observed in 
the laboratory and in space, see e.g. Refs. [28, 291 [loll . 

fall into this category. Q-branch transitions with hKa = and 
transitions with AA'^ - +2 also have considerable intensities 
because of k being very different from -1. The former type of 
transitions was detected in space, the latter both in the labora- 
tory list] and in space. 

The line survey of CW Leo was carried out with the HIFI 
instrument ifssll on board of the Herschel satellite [34] May 
11-15 2010 employing all six bands. Some details on the 
observations and on the data reduction have been given ear- 
lier (iH, and additional information will be provided in a sub- 
sequent manuscript on the whole line survey. The lower reso- 
lution mode provides a point spacing of about 1 . 1 MHz which 
is well sufficient even at low frequencies. A fair fraction of the 
moderately strong to weak lines in the four lowest bands from 
480 to 1120 GHz could be assigned to SiCa. No lines were 
observed with sufficient signal-to-noise ratio in bands 5 to 7 
(1120-1250 GHz and 1410-1910 GHz). The center frequen- 
cies were determined from a fit to the observed line profile using 
the CLASS program of the GILDAS packageH The rotational 
temperature of SiC2 was determined as ~204 K; the Boltzmann 
peak at this temperature is near 500 GHz. 261 distinct features 
have been identified in the HIFI spectra corresponding to 319 
rotational transitions with 19 < 7 < 53 and Kg < 16. All transi- 
tions were /? -branch transitions with AKg - 0. The accuracies 
of the frequencies range from 1.5 MHz for relatively strong, 
isolated lines at low frequencies to 20 MHz for lines with low 
signal-to-noise ratio at high frequencies. Selected SiC2 transi- 
tions from the HIFI molecular line survey are shown in Fig.[T] 

Additional spectral data of IRC +10216 in the radio domain 
(90-358 GHz) were obtained with the IRA^l 30 m telescope. 



located at Pico Veleta, Granada (Spain). Lines in the A 2 mm 
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wavelength domain have been previously reported il37[l . while 
those lying at /I 3 mm and 1 mm will be presented in detail in 
forthcoming papers. The center frequencies were determined 
in the same way as used for the lines observed with Herschel. 
A total of 59 SiC2 lines were clearly observed in emission with 
a spectral resolution of 1-1.25 MHz below 260 GHz and of 
2 MHz above this frequency. The 59 lines correspond to 75 
rotational transitions with Ka < 12. Most of them are R- 
branch transitions with AKg = 0, but there were also five Q- 
branch transitions (AJ = 0) with =2-0 and two wiffi 
Ka =2-2. Although the observed lines are relatively wide, 
around 29 kms ' in equivalent radial velocity, their edges are 
very sharp with a small broadening at the base due to micro 
turbulence velocities of 1.5 kms"'. Hence, lines with a high 
signal-to-noise ratio can be fitted with an accuracy better than 
1 km s"' . Consequently, the experimental accuracies on the fre- 
quencies derived towards IRC +10216 depend on the line in- 
tensity, being around 1 MHz for the weakest lines and around 
0.1 MHz for the strongest ones. 

The transition frequencies from the HIFI and IRAM 30 m 
observations are given in the supplementary material togeffier 
with those from laboratory experiments with quantum numbers, 
uncertainties, and residuals o-c between observed frequencies 
and those calculated from the final fit. 



3. Spectroscopic analysis 

The choice of spectroscopic parameters is often not unique if 
one has to reproduce a large set of transition frequencies, in par- 
ticular in the case of floppy molecules, such as SiC2. One way 
of getting a comparatively small and possibly even unique data 
set is to examine carefully which parameter improves the qual- 
ity of the fit the most, as judged by the rms error, also known 
as reduced x^, and if one includes and keeps only parameters 
determined with a certain level of significance. However, this 
procedure may run into difficulties if parameters of a given or- 
der have rather similar magnitudes or if certain parameters are 
strongly correlated. Moreover, a particularly small parameter 
set or one yielding the smallest rms error does not always yield 
the most reliable predictions, in particular if the amount of tran- 
sition frequencies is rather small. 

As indicated in section [l] the SiC2 laboratory data from 
Ref. were reproduced there to only four times the reported 
uncertainties employing Watson's A-reduction of the rotational 
Hamiltonian with a complete set of parameters up to sixth or- 
der As will be shown later, the quality of the fit is a result of the 
particular choice of parameters, not a result of the choice of the 
reduction. In that fit, the residuals were distributed unevenly, 
and the /iT-ordering of the energy levels was incorrect starting at 
J -26. Refs. [15..,, 17] used a dedicated internal rotation Hamil- 
tonian with 16 parameters and Watson's 5 -reduction of the ro- 
tational Hamiltonian with 17 varied parameters, respectively, to 
reproduce these data within about 1.5 times the reported un- 
certainties. This suggests that the laboratory data ||9[[ were es- 
timated too optimistically. Here, as in previous fits 1 1/7, 36ll . 
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Figure 1: Sections of the molecular line survey of CW Leo displaying selected SiC2 transitions. Their quantum numbers J'^, ^, — J"^,, ^„ are given in red (grey). 
The formulae and quantum numbers of other lines are given in blue (black). 
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the uncertainties from that work were muhipHed with 1 .5 and 
rounded up. 

Starting initially with the parameter set from Ref. |17], the 
complete list of transition frequencies from HIFI observations, 
from laboratory spectroscopy, and from selected additional as- 
tronomical observations could be reproduced well with the S- 
reduction of the rotational Hamiltonian. However, these fits 
required around 24 varied parameters compared with 17 and 
18, respectively, in previous fits [17, Therefore, attempts 
were also made to fit the data employing the A-reduction. In 
one fit, the procedure described above was followed; a fit with 
only 18 spectroscopic parameters achieved an rms error of 0.77. 
However, comparatively many ofF-diagonal distortion parame- 
ters, such as ipK, Ik, and pK, were used in the fit. Moreover, the 
partition function values, calculated by including energy level 
up to y = 90 and Kg - 40 for convergence at 300 K, diff'ered 
remarkably little between 10 K and 300 K. This was caused 
by decreasing energies, eventually turning negative, at low Ka. 
Since such energies are unphysical, the fit was therefore dis- 
carded. 

In a subsequent fit, off-diagonal distortion parameters were 
avoided as far as possible. The fit achieved an overall rms er- 
ror of 0.73 with 20 parameters. Three highly /T-dependent off"- 
diagonal distortion parameters and 0^, used in the previous fit, 
were replaced by I jk and some additional diagonal distortion 
parameters. These parameters yielded partition function values 
similar to previous works flT*, '36'] . The laboratory data from 
Ref. ist] as well as astronomical observations up to 360 GHz 
were reproduced with partial rms errors marginally smaller than 
1 .0, while the transition frequencies obtained from HIFI obser- 
vations had a partial rms error of 0.63, suggesting that, overall, 
they may have been judged slightly conservatively. The param- 
eter (pK, determined insignificantly in the previous laboratory 
spectroscopic work [9] and judged to be rather large in a very 
recent laboratory study of ^^SiCa and ^*'SiC2 , was still not 
determined with significance in later fits and was therefore not 
included in the final fit. The resulting spectroscopic parame- 
ters are given in Table [T] The change in reduction of the ro- 
tational Hamiltonian prompted an attempt to fit the laboratory 
data 1 8, 9] with an A reduction parameter set. All transition fre- 
quencies from astronomical observations were eliminated from 
the line list, and all parameters with uncertainties not much 
smaller than the value were eliminated up to before the point 
at which the rms error increased substantially. The final fit con- 
sisted of 15 parameters which reproduced the data as well as the 
fits described in Refs. [15, 17], but actually with one and two 
varied parameters, respectively, less. The spectroscopic param- 
eters from that fit have also been given in Table [T] 

4. Discussion 

The transitions of SiCi, whose frequencies have been mea- 
sured in the laboratory or in space, are almost all 7J-branch tran- 
sitions, and most of them have lS.Ka - 0. Spectroscopic param- 
eters obtained from such a data set permit reasonable extrap- 
olation up to around twice the highest frequencies in the line 
list even under favorable conditions. Favorable conditions are 



Table 1: Spectroscopic parameters" (MHz) for SiCo obtained in the present 
investigation. 



Parameter 


All data 


Lab. data only'^ 


A - (B + C)/2 


40673.821(37) 


40674.109(60) 


(B + C)/2 


11800.14670 (66) 


11800.14722 (89) 


(B - C)/4 


679.28139 (53) 


679.28365 (87) 




-1.2841(89) 


-1.2092(153) 




1.538195 (69) 


1.538100(86) 


Aj X 10^ 


13.1962(28) 


13.2188 (38) 


Sk X 10' 


869.88(20) 


870.38 (29) 


Sj X 10^ 


2.41187(170) 


2.42028 (374) 


<1>KJ X 10« 


381.0(33) 


426.9(102) 


<1>JK X 10* 


-48.14(81) 


-61.19(308) 


0j X 10' 


-84.9(36) 




IpK X 10^ 


1.084(16) 


0.824(59) 


^JK X 10^ 


-33.51 (43) 


-29.96(113) 


Lkkj X 10' 


319.6(225) 


135.3(207) 


LjK X iC 


-148.4(43) 


-92.9(69) 


LjjK X 109 


-1.43 (31) 


-5.94 (72) 


IjK X 10' 


-1.575(153) 




Pkkj X 10" 


-1.179123) 




Pkj X 1012 


426.3(271) 




PjK X 1012 


-49.50(227) 





" Watson's A-reduction was used in the representation /' . Numbers in paren- 
theses are 1 cr uncertainties in units of the least significant figures. 
Refs. Si 

a good quantum number coverage (several transition frequen- 
cies for each spectroscopic parameter with a sufliciently repre- 
sentative range of J and Kg quantum numbers) and a Hamil- 
tonian that converges sufficiently fast. A predicted transition 
frequency is deemed to be reasonable in this context if the mea- 
sured line is found within three to five times the predicted un- 
certainties. The SiCi molecule is comparatively floppy. More- 
over, the laboratory data set for the main isotopolog of SiCa is 
rather sparse and extends to 370 GHz. Nevertheless, the high- 
est frequency transition in Ref. li36ll . which report a preliminary 
analysis of HIFI data, is 30o,3o - 29(),29 at 636346 MHz, and it 
shows among those lines the largest deviation from the predic- 
tions based on Ref. [i 17] with almost 10 MHz. However, this de- 
viation corresponds to only 5.5 times the predicted uncertainty. 
Moreover, the frequency from the final analysis has been cor- 
rected to a slightly lower frequency, in better agreement with 
the initial predictions. 

Larger deviations from the predictions were observed around 
1 THz, but these had to be expected since the predicted tran- 
sition frequencies have uncertainties of several tens to almost 
100 MHz. In fact, predictions from Refs. [17, 361 are quite rea- 
sonable for the most part even if the deviations are large enough 
to be easily recognizable in the astronomical spectra. 

Fits employing Watson's A-reduction of the rotational 
Hamiltonian reproduce the experimental transition frequencies 
with fewer parameters than those using the 5 -reduction and 
should thus be preferred as long as the energy ordering is 
correct up to sufficiently high quantum numbers. It may be 
interesting to compare the quality of the predictions of the 
transitions detected with HIFI from either laboratory fit from 
Refs. linn or in Table [1] Predictions based on 5 -reduction pa- 
rameters f\r\ are quite good for transitions having higher val- 
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ues of Ka (6 - 14), but deteriorate for lower Kg values. In con- 
trast, predictions based on A-reduction parameters (Table[T]un- 
der heading "Lab. data only") are quite good for low-^Tf, tran- 
sition, but deteriorate for higher Ka values, which is closer to 
common expectations. Both fits as well as both predictions, 
as well as additional fits are available in the archive section of 
the CDMSHI Overall, however, the quality of the predictions is 
quite similar and do not provide a reasoning which reduction 
should be preferred. A very recent laboratory study of ^^SiCa 
and ^"SiCa 11 111 also obtained fits of similar quality for the A- 
and the S -reduction, but only if one more parameter was used 
in the 5 -reduction and preferred thus the A-reduction. 

The SiCa spectroscopic parameters in Table [1] from labora- 
tory data only and from all data agree quite well as far as they 
have been determined, in particular for the lower order param- 
eters, rotational and quartic distortion constants. The higher 
order parameters show larger deviations, possibly reflecting the 
floppy nature of the molecule. However, the deviations are in 
all instances, at most, barely significant, meaning larger than 
three times the initial uncertainties. 

The measured transition frequencies have been reproduced 
well with parameter sets of conventional Watson-type Hamilto- 
nians. The molecule is rather floppy already in its ground vibra- 
tional state such that a fairly large number of parameters were 
used. Additionally, the magnitudes of the parameters decrease 
slowly such that the contributions of the higher order parame- 
ters are fairly large at energy levels having high values of J or 
Ka- Extrapolations to slightly higher quantum numbers may be 
reasonable, but should be viewed with caution. 

The dedicated internal rotation Hamiltonian OlSll did not 
achieve a considerably more compact representation of the 
measured transition frequencies known at that time. Without 
any available predictions for transition frequencies or rotational 
energies it seems questionable if that model is more appropri- 
ate than a conventional Watson-type Hamiltonian, even more so 
as the predicted barrier to linearity was rather far from values 
from ab initio calculations [|7t] or from the SRB extrapolation 
from highly vibrationally excited levels [16]. Moreover, one 
may question if an internal rotation model is appropriate at all. 
On the other hand, the semirigid bender approach [ 16] may be 
a viable alternative, but only if it will achieve reproduction of 
the measured transition frequencies to within estimated uncer- 
tainties with a parameter set that is not much larger than those 
used previously. 

Additional laboratory data for the main isotopolog of SiC2 
will be of great use in particular if these connect levels with 
Ka >2 via transitions with AKa > 2. But even Q-branch tran- 
sitions with AKa - may be very useful as they sample higher 
J levels at comparatively low frequencies. Such higher excited 
transitions may well be observable with existing interferometric 
facilities such as PdBl, ATCA, or SMA, and even more so with 
upcoming instruments such as ALMA, NOEMA, or possibly 
EVLA. 



5. Conclusions 

Astronomical observations with Herschel have been used to 
improve spectroscopic parameters of SiC2 in its ground vi- 
brational state greatly. The transition frequencies obtained 
from laboratory data as well as from astronomical observa- 
tions have been reproduced to within the estimated uncertain- 
ties employing a conventional Watson-type Hamiltonian. These 
data should be useful to test alternative models to describe 
the rotational or even rovibrational energy levels of this floppy 
molecule. Predictions of the rotational spectrum as well as line, 
parameter and other auxiliary files, both from present fits as 
well as previous ones, will be available in the CDMS 1261. 271. 
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Appendix A. Supplementary Material 

Supplementary data for this article are available on Sci- 
enceDirect (www.sciencedirect.com) and as part of the 
Ohio State University Molecular Spectroscopy Archives 
(http://library.osu.edu/sites/msa/jmsa_hp.htm). Supplementary 
data associated with this article can be found, in the online ver- 
sion, at doi: . 
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